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ABSTRACT 

We present an investigation of the properties and environments of bright extremely 
red objects (EROs) found in the fields of the quasars TXS 0145+386 and 4C 15.55, both 
at z ~ 1.4. There is marginal evidence from Chandra ACIS imaging for hot cluster gas 
with a luminosity of a few 10^^ ergs s~^ in the field of 4C 15.55. The TXS 0145+386 
field has an upper limit at a similar value, but it also clearly shows an over density of 
faint galaxies. None of the EROs are detected as X-ray sources. For two of the EROs 
that have spectral-energy distributions and rest-frame near-UV spectra that show that 
they are strongly dominated by old stellar populations, we determine radial-surface- 
brightness profiles from adaptive-optics images. Both of these galaxies are best fit by 
profiles close to exponentials, plus a compact nucleus comprising ~ 30% of the total 
light in one case and 8% in the other. Neither is well fit by an r^/^-law profile. This 
apparent evidence for the formation of massive ~ 2 x 10^^ disks of old stars in the early 
universe indicates that at least some galaxies formed essentially monolithically, with 
high star-formation rates sustained over a few 10* years, and without the aid of major 
mergers. 
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Subject headings: galaxies: high-redshift — galaxies: formation — galaxies: evolution 

1. Introduction 

Galaxies classified as extremely red objects (EROs; e.g., Daddi et al. 2000; Chapman et al. 2000; 
Stockton 2001; Cimatti et al. 2002; Saracco et al. 2004; Yan, Thompson, k Soifcr 2004; Daddi et 
al. 2004; earlier references can be found in these papers) are usually defined to have R — K greater 
than a value somewhere in a range of 5-6. They have long been known to include at least two 
disparate classes. Some are highly reddened high-redshift star-forming galaxies (e.g., Dey et al. 
1999; Cimatti et al. 2002; Yan et al. 2004b). The remainder are almost exclusively galaxies with 
2; > 1 whose stellar populations arc already ancient and have little admixture from more recent star 
formation (e.g., Stockton, Kellogg, & Ridgway 1995; Dunlop et al. 1996; Spinrad et al. 1997; Soifer 
et al. 1999; Liu et al. 2000; Stockton 2001; Stockton, Canalizo, & Maihara 2004; Chen & Marzke 
2004; Cimatti et al. 2004; Fu, Stockton, & Liu 2005). For convenience, we shall henceforth simply 
refer to this latter class as "old galaxies" (OGs; however, note that this term refers to the age of the 
stellar population, not necessarily the age of assembly of the galaxy itself). These OGs represent 
the earliest major episodes of star formation in the universe, and they seem to have formed the 
overwhelming majority of their stars over a very short period. Accordingly, they potentially offer 
us an avenue to investigate the formation and early evolution of the very first massive galaxies to 
appear in the universe. 

There have been a number of discussions of the mix of morphologies in various samples of 
EROs. It is commonly assumed that morphological types are surrogates for the two classes of 
EROs mentioned above; i.e., that spheroidals as a group are more-or-less congruent with galaxies 
dominated by old stellar populations and that disk-like galaxies and those showing obvious merging 
activity arc dusty star-forming galaxies (Moriando, Cimatti, & Daddi 2000; Stiavelli & Treu 2001; 
Smith et al. 2002; Sawicki ct al. 2005). Nevertheless, both Gilbank et al. (2003) and Yan, Thompson, 
& Soifer (2004) mention the possibility of truly passive disk-dominated galaxies. lye et al. (2003) 
seem to have found such a galaxy at z ^ 1.5, and Stockton, Canalizo, &: Maihara (2004) and Fu, 
Stockton, & Liu (2005) have shown examples of galaxies that seem to be essentially pure disks of 
old stars at 2; = 2.48 and z = 1.34, respectively. 

So far, the morphologies of only a few OGs at redshifts significantly greater than 1 have been 
studied in any detail (e.g.. Bunker et al. 2003; Stockton, Canalizo, &; Maihara 2004; Yan, Thompson, 
& Soifer 2004; Cimatti et al. 2004; Fu, Stockton, & Liu 2005; Daddi et al. 2005). Here we present 
an investigation of the morphologies and environments of 3 OGs (and 1 dusty star- forming galaxy) 
found in the fields of the quasars TXS 0145+386 (z = 1.442) and 4C 15.55 {z = 1.406). We assume 
a fiat cosmology with Hq = 70 and O^, = 0.3. 
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2. Observations and Data Reduction 

2.1. Optical and Near-IR Observations 

The two fields discussed in this paper are part of a larger sample identified from a survey of 
the fields of (radio-loud) quasars with 1.4 < z < 1.6. We have chosen to observe quasar fields in 
this redshift range for three reasons. First, there is the simple and practical observational matter 
that, for our near-IR observations, it is necessary to have one or more objects in each field that 
will be detected in single short exposures and that can be used to center each dithered frame that 
is coadded to make the final image. Centering the field on a quasar ensures that at least one such 
object will be present, even when our field size is relatively small. Second, radio sources at high 
redshifts appear to mark regions that include those of high overdensity (e.g.. Best 2000; Best et al. 
2003; Barr et al. 2003; Rottgering et al. 2004), which are precisely the regions in which most models 
expect processes of galaxy formation and evolution to have proceeded most rapidly. Significantly, 
compelling results from the 2dF survey show that radio quiet QSOs are much less strongly clustered 
than are quasars (Croom et al. 2001). Finally, and most importantly, looking for companions to 
radio sources at a specific redshift allows us to choose redshifts for which the photometric diagnostics 
from standard broadband filters give the cleanest separation between old galaxy populations and 
highly reddened star-forming galaxies or other possible contaminants. Specifically, for z ~ 1.5, the 
4000 A break falls just shortward of the J band. At this redshift, an R— J or I— J vs. J—K plot 
provides excellent discrimination between OGs and most other objects with similar R—K values 
(Pozzetti & Mannucci 2000; Stockton 2001). Cimatti et al. (2003) claim that some star-forming 
galaxies may have colors similar to those of OGs; however, their Fig. 5, on which this claim is 
based, seems to depend strictly on morphological classes judged by eye from HST ACS images, 
rather than on spectroscopic evidence for star formation. This approach assumes that the generally 
tight correlation between morphology and star-formation rate that applies at the present epoch 
continues to hold at high redshifts. One of the purposes of the present paper is to discuss evidence 
that this assumption may not always be valid. 

We have used a photometric sieve process to isolate candidates. The initial observations were 
carried out at K' with NSFCam on the NASA Infrared Telescope Facility (IRTF). If any objects 
were found within a 30" square region surrounding the quasar (limited by the NSFCam field and 
the dithering amplitude) and with 17 < K' < 19 (we use the Vega magnitude scale throughout this 
paper), we next obtained J-band imaging, li J—K' k.2 (as expected for old stellar populations in 
this redshift range), we then obtained R (and sometimes I and/or z') imaging, generally with the 
Echellette Spectrograph and Imager (ESI) on the Keck II telescope, but sometimes with other CCD 
systems and telescopes. Since we were searching for systems of old stars with essentially no recent 
star formation, we required R — K' > 6; this is, of course, by itself, a necessary but not sufficient 
condition at redshifts z ^ 1.5 for the absence of young stars, because of the possible effects of 
reddening. However, together with the constraint that J—K' « 2, this R—K' color requires a sharp 
inflection in the spectrum consistent with the 4000 A break for old stellar populations at z ~ 1.5. 
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We show traces for an unreddened, passively evolving galaxy and for various models of heavily 
reddened galaxies as functions of redshift in the J — K' — R—J diagram in Fig. 1, along with the 
positions in this diagram of the 4 objects discussed in this paper. The passively evolving galaxy, 
deliberately chosen to be an extreme case, is assumed to have formed all of its stars 5 x 10^ years 
after the Big Bang. Reddened starbursts are represented by constant-star-formation-rate models, 
with Calzetti et al. (2000) and Ferrara et al. (1999) reddening adjusted to give R — K' = 6.0 at 
z = 1.4. The reddening curve of Ferrara et al. (1999) used here is that of a piuc disk with embedded 
dust and an inclination of 70° to the line of sight. Finally, we also show a model for continuous, 
exponentially decreasing star formation, with an e-folding time of 5 Gyr and an age of 12 Gyr, again 
subjected to Calzetti et al. (2000) reddening to give R—K' = 6.0 at z = 1.4. The age of this model 
is of course unrealistic at the redshifts we are considering, but it serves as a proxy for a range of 
models for which there is a substantial old stellar component coupled with a relatively modest star- 
formation rate and high reddening. We will present observations of the larger sample elsewhere. 
We restrict ourselves here to the TXS 0145+386 and 4C 15.55 fields because we have obtained more 
complete observations for these fields than for any of the others, including spectroscopy of one 
OG in each field with the Low- Resolution Imaging Spectrograph (LRIS; Oke et al. 1995) on Keck 
I and/or ESI on Keck II. The spectroscopy with LRIS always used a 1" slit or slitmask with a 
400 line/mm grating blazed at 8500 A, and the ESI spectroscopy always used a 1''25 slit and the 
echellette mode. Other follow-up observations on the TXS 0145+386 and 4C 15.55 fields include 
deep imaging with the University of Hawaii 88-inch Telescope (UH88) and the United Kingdom 
Infrared Telescope (UKIRT), and adaptive optics (AO) imaging with the Canada- France-Hawaii 
Telescope (CFHT) and the Subaru Telescope (Kobayashi et al. 2000). For the galaxy morphology 
analysis from AO images, we have obtained PSF information from the closest sufficiently bright 
star in the field (12" away for TXS 0145+386 ERl and 9" away for 4C 15.55 ER2). By mapping out 
the PSF variation with position from other stars in the field and from other fields obtained on the 
same nights, we confirm that there should be no significant change in the PSF over these distances 
on the nights of our AO observations. 

A listing of all of the optical/IR observations relevant to this paper is given in Table 1. 

The optical and IR imaging observations were processed through our standard IRAF reduction 
pipelines. The spectroscopic data needed more detailed and individual attention, both in the 
planning of the observations and the reduction, because of the faintness of the objects. For the 
longslit observations, both with LRIS and with ESI, we arranged a slit position angle that would 
allow a brighter object to be placed near one end of the slit while the actual target was more-or- 
less centered. After the individual spectra were flat-fielded, distortion corrected, and wavelength 
calibrated, the spectrum of the bright object could be used to register the spectra for coadding 
before extracting the spectrum of the target object. For the ESI observations, the bright object 
was also necessary to confirm that the echellette orders were properly normalized to each other in 
the overlap regions before they were combined to make a continuous spectrum. 



-5- 



Table 1. Optical and IR Observations of the TXS 0145+386 and 4C 15.55 Fields 



Date (UT) 


Telescope 


Iiistruiiieiit 


Target 


Filter 


Exposure 


Pixel Scale 


Comments 


1998 Aug 17 


IRTF 


NSFCam 


4C 15.55 Field 


K> 


540 


s 


0^30 




1998 Aug 17 


IRTF 


NSFCam 


TXS 0145+386 Field 


K' 


540 


s 


0'.'30 




1998 Aug 18 


IRTF 


NSFCam 


4C 15.55 Field 


J 


540 


s 


O'.'SO 




1998 Aug 18 


IRTF 


NSFCam 


TXS 0145+386 Field 


J 


540 


s 


0^30 




1998 Aug 19 


IRTF 


NSFCam 


4C 15.55 Field 


K' 


540 


s 


0^'30 




1998 Aug 19 


IRTF 


NSFCam 


4C 15.55 Field 


J 


1620 


s 


0^'30 




1998 Aug 19 


IRTF 


NSFCam 


TXS 0145+386 Field 


K' 


1080 


s 


0^'30 




1998 Aug 19 


IRTF 


NSFCam 


TXS 0145+386 Field 


J 


1080 


s 


0^'30 




1998 Sep 05 


UH88 


QUIRC 


TXS 0145+386 Field 


J 


3000 


s 


0'.'18 




1998 Oct 28 


UH88 


Tek2048 


TXS 0145+386 Field 


Rc 


12000 


s 


0^'22 




1998 Oct 29 


UH88 


Tek2048 


TXS 0145+386 Field 


Ic 


14400 


s 


(y.'22 




1999 Aug 24 


CFHT 


Pueo/KIR 


TXS 0145+386 Field 


K' 


6000 


s 


0'.'035 




1999 Oct 02 


Keck I 


LRIS 


TXS 0145+386 Field 


Z 


16200 


s 


0'.'22 




1999 Nov 03 


Keck I 


LRIS 


TXS 0145+386 ERl 


Spec 


16800 


s 


0'.'22 


1" slit 


2000 Jun 05 


Keck II 


ESI 


4C 15.55 ER2 


Spec 


18000 


s 


0'.'15 


1'.'25 slit, echellette mode 


2000 Jun 05 


Keck II 


ESI 


4C 15.55 ERl 


Spec 


3600 


s 


0^'15 


1''25 slit, echellette mode 


2000 Oct 03 


Keck II 


ESI 


TXS 0145+386 ERl 


Spec 


21600 


s 




1"25 slit, echellette mode 


2000 Oct 04 


Keck II 


ESI 


TXS 0145+386 ERl 


Spec 


19120 


s 


(y.'ib 


1''25 slit, echellette mode 


2001 Aug 24 


Keck I 


LRIS 


TXS 0145+386 ERl 


Spec 


10800 


s 


0^22 


slitmask 


2002 May 31 


Subaru 


CISCO 


4C 15.55 Field 


K' 


1920 


s 


O'.'ll 




2002 May 31 


Subaru 


CISCO 


4C 15.55 Field 


J 


960 


s 


O'.'ll 




2002 Aug 07 


Keck II 


ESI 


4C 15.55 Field 


Rc 


2520 


s 


0'.'15 


Ellis R filter 


2002 Aug 17 


Subaru 


IRCS/AO 


4C 15.55 Field 


K' 


2160 


s 


0^'058 




2003 Dec 05 


UKIRT 


UFTI 


TXS 0145+386 Field 


K 


8400 


s 


(y.'o9o 
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Fig. 1. — Positions of model tracks and the objects discussed in this paper in the J — K' vs. R—J 
plot. Each of the model tracks ranges from z = 0.4 to z = 2.0 (bottom to top or left to right); 
increments of 0.1 in Az are indicated by the black dots, with the z = lA point circled. The 
red trace is that of a passively evolving stellar population for which all of the stars were formed 
instantaneously when the universe was 5 x years old (corresponding to z = 9.5). The blue 
trace shows a continuous star-formation model with an exponentially decreasing star-formation 
rate having an age of 12 Gyr and an e-folding time of 5 Gyr, with sufficient Calzetti et al. (2000) 
reddening to produce R — K' = 6.0 at z = 1.4. The green and magenta traces show models for 
starbursts, with a constant star- formation rate and Calzetti et al. (2000) and Ferrara et al. (1999) 
reddening, respectively, both also normalized to give R — K' = 6.0 at 2; = 1.4. See the text for 
details. 

2.2. X-Ray Observations 

The purpose of the Chandra X-ray observations was to search for evidence for cluster gas 
surrounding the two quasars. The TXS 0145-1-386 and 4C 15.55 fields were observed with ACIS-I on 
Chandra on January 4, 2002 UT for 63 ksec, and on April 13-14, 2002 UT for 52 ksec, respectively. 
ACIS-I consists of four CCDs with a total field of view of 16' x 16', covering an energy range from 
0.3 to 10 keV. At full resolution, the plate scale is 0'.'5 per pixel. The images were obtained in the 
VFAINT mode, which potentially offers better rejection of spurious events and therefore provides a 
lower background, at the possible expense of exacerbating pileup problems near bright sources. We 
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analyzed both the pipehne-reduced images, which used the standard algorithm for grading events, 
and images we reprocessed using the VFAINT algorithm. 

Using the Chandra Interactive Analysis of Observations (CIAO) software^, we produced expo- 
sure maps for each observation which wc then used to flat-field the images and convert to flux units. 
In order to simplify the calculation of the exposure map, we assumed a monoenergetic distribution 
of source photons for both observations. Ideally one would prefer to define a set of spectral weights 
for the image in order to more accurately determine the exposure map, which is a complicated 
function of energy. However, given that each of the sources are strongly peaked at energies of 
~ 1.1-1.3 keV, and that exposure maps produced at slightly higher and lower energies produced 
negligible difference in the final fluxed images, we are confident that the error introduced by using 
the monoenergetic approximation is small. 

For each observation, source and background regions were chosen from the final fluxed image 
and then applied to the raw counts image. Several background regions were chosen for statistical 
robustness, and the dispersion in the background regions was found to be approximately 17% 
for TXS 0145+386 and 14% for 4C 15.55; these percentages are similar to the calculated Poisson 
errors. Prom the counts image, we calculated the Poisson upper and lower confidence limits for the 
background and source regions using the values given by Gehrels (1986) and Ebeling (2003). 

We computed model PSPs for each source using the Chanda Ray Tracer (ChaRT) tool, which 
contains detailed information about the High Resolution Mirror Assembly (HRMA) on Chandra, 
and allows for highly accurate determinations of the PSF over a range of energies and off-axis 
positions. Using CIAO's modeling and fitting program, Sherpa, we were able to model the spectrum 
of each source with a 1-D power law combined with photoelectric absorption. This spectrum was 
passed to ChaRT in order to determine the ray traces, which were then projected onto the ACIS-I 
detector plane, and converted to a PITS image using MARX. Both TXS 0145+386 and 4C 15.55 are 
bright sources and suffered from photon pileup; however, MARX allowed us to simulate the PSP 
both with and without pileup. We find that while pileup affects the innermost pixels in the PSP, 
the outer wings of the PSF (> 5") are essentially unaffected, even with the VFAINT processing, so 
that an accurate subtraction of the contribution from the PSP can still be obtained in the region 
relevant to our search for cluster emission. 

3. Results 

3.1. The TXS 0145+386 Field 

Images in several bandpasses of the TXS 0145+386 field are shown in Fig. 2. Two galaxies 
that meet our selection criteria are found in this field; these are labelled "ERl" and "ER2" in 



^http: / /cxc.harvard.edu / ciao/ 
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the figure. The spectral-energy distributions (SEDs) for these galaxies are shown in Fig. 3, along 

with a reference spectrum for a 3-Gyr-old solar-metallicity instantaneous-burst model (Bruzual &: 
Chariot 2003). There appears to be an excess of faint objects in Fig. 2 surrounding the quasar 




Fig. 2.— Images of the TXS 0145+386 field. The quasar TXS 0145+386 itself is indicated by the 
"Q". The two EROs are labelled "ERl" and "ER2" and are marked by the short lines in the K 
image. 

and ERl. This impression is confirmed by Fig. 4, which shows a 4f3-square region centered on 
TXS 0145+386, with positions of all galaxies in a magnitude interval 24 < z < 25.5 indicated. We 
have attempted to obtain spectroscopy of TXS 0145+386 ERl. Figure 5 shows that the spectrum 
is generally consistent with that expected for an old stellar population at the approximate redshift 
of the quasar. However, since the S/N ratio for this spectrum is not very high and because there 
are no strong absorption features in the available spectral range, it is also probably consistent with 
the average "dusty star-forming" spectrum of Cimatti et al. (2002, see their Fig. 2). The presence 
of significant star formation would seem, at least at first sight, to be corroborated by the presence 
of [O II] A3727 emission. However, spectra obtained at higher resolution with ESI on Keck II show 
that the emission has an intrinsic velocity width of 900 km s~^, considerably higher than would 
be expected simply from ordered rotation or other gravitationally induced motions of star-forming 
regions in even the most massive galaxies. Such velocities could indicate the presence of wind. 
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Fig. 3. — Photometry of two EROs found in the TXS 0145+386 field. The photometric data are 
shown in black, with the vertical bars indicating Icr uncertainties and the horizontal bars indicating 
filter half-transmission widths. The green curves show a Bruzual & Chariot (2003) 3-Gyr-old 
instantaneous burst model for reference. 

driven either by a hidden starburst or by a hidden AGN. On the other hand, as we describe in 
more detail below, there is no evidence for a point X-ray source that would indicate the presence of 
an AGN at this location in our deep Chandra ACIS images. The absence of obvious Mg II or [Ne V] 
emission also favors stellar excitation over an AGN origin. If the emission is due to recent star 
formation, the [O II] luminosity indicates a current star-formation rate (uncorrected for internal 
extinction) of ~ 2 Mq yr~^, using the conversion given for high-redshift galaxies by Kewley, Geller, 
& Jansen (2004, eq. 19). This rate is very similar to the rate found by Cimatti et al. (2002) for 
their sample of dusty star- forming EROs. 

In sum, the [O II] emission seems most reasonably interpreted as indicating a moderate degree 
of star formation (depending, of course, on the amount by which the [O II] emission is affected by 
internal extinction); the large velocity width remains a puzzle. However, the SED (Fig. 3), with 
its strong break at redshifted 4000 A, clearly matches that of an old stellar population. We know 
of no way to reproduce a similar SED with any reasonable reddening of a stellar population whose 
rest-frame UV-optical luminosity is dominated by young stars. It seems safe to conclude that the 
observed-frame K'-hand morphology of TXS 0145+386 ERl is overwhelmingly dominated by light 




Fig. 4. — Evidence for clustering of faint galaxies in the vicinity of TXS 0145+386. The green dots 
show the distribution of galaxies with 24 < z < 25.5 in a 4'.3 x 4'.3 field around TXS 0145+386. 
Black areas are regions obscured by bright stars or galaxies. The two EROs are indicated by red 
crosses, and the quasar by a blue cross. 
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Fig. 5. — The spectrum of TXS 0145+386 ERl (black line), obtained with LRIS. The red trace 
shows the Bruzual &: Chariot (2003) 3-Gyr instantaneous-burst model. Note the [O II] A3727 
emission in the galaxy spectrum, which gives a redshift z = 1.4533 for the galaxy. 
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from old stars. 

There is a mag 12 star about 1' south of TXS 0145+386 ERl. Although this offset is rather 
larger than one would like, we were able to obtain AO imaging with PUEO on CFHT on a night 
during which the angular coherence was exceptionally good, as indicated by the round PSFs of the 
quasar and the nearly stellar galaxy between the quasar and ERl. The K' AO image is shown 
in Fig. 6. TXS 0145+386 ERl itself shows an elliptical image with a high degree of symmetry, 
surrounded by some irregular structure at very low surface brightness. Determining the radial 




Fig. 6.— AO imaging of TXS 0145+386 ERl with PUEO and KIR on CFHT. In the left panel, the 

quasar is at the bottom, ERl at the top. The object near the center is an extremely compact galaxy 
shown by our spectroscopy to be at z = 0.8. The right panels show, from top to bottom, two differ- 
ent contrast stretches for each of three degrees of smoothing of enlarged images of TXS 0145+386 
ERl. Although there are clearly problems with uneven background levels in the higher contrast 
images, there seems to be some evidence for irregularities in the low-surface-brightness outer parts 
of ERl. 

surface-brightness profile for TXS 0145+386 ERl from the AO image is complicated by the fact 
that there is an unresolved core present, which is offset from the centroid of the underlying galaxy 
by ~ 0'.'12 (= 1 kpc). Because both the relative centroids and the ampUtudes of the core and 

galaxy components vary slightly in different models, wc have carried out independent modeling of 
the galaxy in terms of r^/^-law, exponential, and Sersic profiles, using GALFIT (Peng et al. 2002). 
In each case, we included a PSF component as well to represent the star-like core, and we allowed 
the centers of the galaxy and PSF components to float freely. Figure 7 shows fits of PSF-convolved 
radial-surface-brightness plots for each of the models to the measured profile, after removal of the 
compact nucleus. The best-fit Scrsic model has a quite low index of 0.34. Even so, the points at 
large radius still tend systematically to fall below the curve. The exponential model is nearly as 
good a fit, and the r^/^-law model is clearly worse, being too high at small and large radii and too 
low at intermediate radii. For the Sersic model, the compact core comprises about 30% of the total 
light in the K' band. 
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r (arcsec) 

Fig. 7. — Seeing-convolved radial surface-brightness plots fitted to the observed radial-surface- 
brightness profile of TXS 0145+386 ERl. The surface-brightness scale to the left is correct for the 
bottom (red) curve and points; the other curves and their corresponding points are displaced by 2 
and 4 magnitudes. The compact nuclear component has been removed from both the models and 
the data. Because the model centers for both the galaxy and the nucleus vary slightly among the 
various models, as does the ellipticity distribution (as a function of radius), we have used these 
model centers and ellipticities to constrain the determination of the observed profiles in each case. 

We do not have AO imaging of TXS 0145+386 ER2; nevertheless, our best K' images show a 
somewhat asymmetric structure suggesting that this galaxy may be an ongoing or recent merger. 
If so, the fact that the SED is clearly that of an old stellar population with no hint of recent star 
formation indicates that the participating galaxies are likely to have been largely free of cold gas. 

For the X-ray statistics, we chose to use an inner radius of 10" and an outer radius of 30" for 
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our annulus around TXS 0145+386. The outer limit was chosen at this scale because of a nearby 
X-ray source that would have confused the statistics had we used a larger radius. This outer radius 
corresponds to a physical scale of 253 kpc at z = 1.442. After subtracting background counts and 
the modeled wings of the PSF in the annulus, we are left with an insignificant 3 (+12, —9) excess 
counts. In order to estimate what this result means in terms of limits on cluster X-ray emission, 
we assume a specific reasonable model: a modified King profile (so-called /3 model) with /3 = 0.7, 
a core radius of 140 kpc, a metallicity of 0.4 solar, and a temperature kT = 6 keV (cf. Mullis et 
al. 2005), and centered on the quasar. Taking all of the flux out to 10 core radii, we find that the 
annulus contains 24.5% of this total, so the extrapolated total residual counts would be 12 (+49, 
—37). We use WebPIMMS^ to calculate the total luminosity corresponding to a 3 a upper limit of 
approximately 150 counts (2.4 x 10~^ counts s~^), assuming a Raymond-Smith model for the SED. 
This limit comes out to 5.7 x 10^^ ergs for a rest-frame energy range from 0.5-8 keV. This limit 
is not a very strong constraint; for example, the cluster recently found by Mullis et al. (2005) at 
a similar redshift has 3 x 10^^ erg s~^ over the 0.5-2 keV range, which, given the difference in 
energy range, is similar to our upper limit. 

There is no significant evidence for X-ray emission from either TXS 0145+386 ERl or TXS 0145+386 
ER2. Formally, for ERl, we have 3 + 3 x 10~^ counts per second, and, for ER2, + 3 x 10^^ counts 
per second. In contrast, for TXS 0145+386 itself, we have 1.75 + 0.05 x 10~^ counts per second. 

3.2. The 4C 15.55 Field 

The field of the quasar 4C 15.55 is shown in Fig. 8. Again, there are two EROs in the field. 
However, our photometry (Fig. 9) indicates that ERl, though quite red {R—K' ~ 6.5), does not 
show the strong break characteristic of an old stellar population (see also Fig. 1); we therefore 
classify it as a likely dusty star-forming galaxy. On the other hand, ER2 in this field clearly meets 
our criterion for an old stellar population. The rather large uncertainty shown for the i?-band 
photometry of ERl in Fig. 9 results from having to correct for the effect of the blue object about 
2" south-east of ERl. Figure 10 shows our spectrum of 4C 15.55 ER2, obtained with ESI on Keck 
II. The S/N of this spectrum is sufficiently high that we can see clear matches between absorption 
features in the spectrum and those in the 3-Gyr Bruzual & Chariot (2003) model. The galaxy has 
a redshift of 1.412, slightly higher than that of the quasar, which is given as 1.410 by Smith et 
al. (1978) and as 1.406 by Wills & Lynds (1978). Our AO imaging of the 4C 15.55 field is shown 
in Fig. 11. In order to have a large enough field to image both EROs in this field, we used the 
0'.'058 pixel scale option. ERl, which we had inferred on photometric grounds to be a reddened 
star-forming galaxy, is clearly an interacting system, while ER2 shows a very regular, symmetric 
profile. 



''http://heasarc.gsfc. nasa. gov /Tools/ wSpimms . html 
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Fig. 8. — Images of the 4C 15.55 field. The quasar 4C 15.55 itself is indicated by the "Q", and 
the two EROs are labelled "ERl" and "ER2". Note the difference in morphology between ERl, 
which the photometry indicates is likely a dusty star- forming galaxy, and ER2, which clearly has 
the characteristics of an old stellar population with little reddening (see Fig. 9). The left and right 
versions of the J and K' images show different contrast levels. 
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Fig. 9. — Photometry of two EROs found in the 4C 15.55 field. The photometric data are shown 
in black, with the vertical bars indicating la uncertainties and the horizontal bars indicating filter 
half-transmission points. The solid green curves show scaled Bruzual Chariot (2003) 3-Gyr-old 
solar-met allicity instantaneous burst models for reference. 
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Fig. 10. — The spectrum of 4C 15.55 ER2 {black line). The red trace shows the Bruzual & Chariot 
(2003) 3-Gyr 0.32 solar-metalhcity instantaneous-burst model. The redshift of 4C 15.55 ER2 is 
1.412. 




Fig. 11. — Images of 4C 15.55 ERl and ER2 obtained with the Subaru AO system and IRCS. 
The left-hand panels show the images smoothed with a Gaussian kernel with a = 0.5 pixel at 
low contrast; the middle and right-hand panels are at successively higher contrasts smoothed with 
a = 1.0 pixel and a = 2.0 pixels, respectively. 
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We have carried out two-dimensional, PSF-convolved modeling of the profile of ER2, using 
GALFIT and the STSDAS ellipse task. As Fig. 12 shows, a pure r^/^-law is a poor fit to the profile. 
An exponential model provides an excellent fit beyond a semi-major axial radius of about 0'.'2 but 
falls below the observed profile at the center. A Sersic profile with n = 2.0 comes reasonably close to 
fitting the points at all radii, although it is slightly above at both large and small radii and slightly 
below in the middle. Our best fit is achieved with a composite profile, comprising a dominant 
exponential with the addition of a small stellar component to the nucleus (this "stellar nucleus" 
could be a r^/^-law bulge or a nuclear disk with a small effective radius; in any case, the contribution 
to the total light from this compact central component is small). For the Chandra ACIS data for 
4C 15.55, the analysis is made difficult by the presence of an X-ray jet associated with the QSO 
(see Fig. 13). This jet likely contributes to the extended emission around the QSO, but, while the 
contribution from the central PSF can be modeled, the jet's full contribution cannot be known 
precisely. We therefore tried to avoid any region where jet emission would confuse the statistics. 
We chose a circular aperture of radius 50" around the QSO (corresponding to '^700 kpc at z = 1.5), 
excluding an elliptical region with a semi-major axis also of 50" and a semi-minor axis of 15" and 
with a position angle of — 68?2, corresponding to the observed jet position angle. We measured 
the flux in background regions using the same aperture and calculated the Poisson errors. After 
accounting for the contribution from the quasar's PSF and subtracting the background, we find an 
excess of 37ibl6 counts around the source. We note that the interpretation of this detection remains 
somewhat ambiguous, given the presence of the jet, but that an overdensity of events is clearly 
visible around the source in the original Chandra image. On the assumption that this overdensity 
can be attributed to an extended cluster component rather than diffuse radiation associated with 
the X-ray jet, we estimate the X-ray luminosity of the cluster using the same assumptions as we 
did for obtaining an upper limit for the TXS 0145-1-386 field. In this case, extrapolating the counts 
in our aperture to 10 core radii gives a total of 93 ± 39. Again using WebPIMMS, with the same 
parameters as before, we find a total luminosity of 4.0 it 1.7 x 10^'^ ergs s~^. This is again similar to 
the luminosity found by Mullis et al. (2005) for their non-AGN-related cluster at the slightly lower 
redshift of 1.39. 

As for TXS 0145-F386 ERl and ER2, 4C 15.55 ER2 itself shows no evidence for point-like X-ray 
emission. 



4. Discussion 

The clustering amplitude of luminous galaxies comprising old stellar populations at high red- 
shifts (Daddi et al. 2000; McCarthy et al. 2001; Brown et al. 2005) is similar to that of early- type 
galaxies at the present epoch. It is also comparable to that of powerful radio sources at similarly 
high redshifts, so it would not be surprising if these two classes were to have a considerable overlap 
in their spatial distribution. Hall et al. (2001) and Wold et al. (2003) have found an excess of EROs 
in radio-loud quasar fields that tends to corroborate such an overlap. Wold et al. (2003) suggested 
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Fig. 12. — The surface-brightness profile of 4C 15.55 ER2, from AO imaging obtained with the 
Subaru AO system and IRCS. The colored curves show fits of various functions to the data, all 
convolved with the point-spread function. The green curve in the upper panel shows an r^/^-law fit, 
normalized to the central data points, but falling below them at intermediate radii and above them 
at large radii. The red curve shows an exponential profile with re = C/.'SG, which is an excellent fit 
beyond 0'.'2, but under-fits the central region. In the lower panel, the green curve shows the best-fit 
Sersic profile, with n = 2.0, and the red curve shows a best-fit composite profile, adding an 8% 
contribution from a stellar or compact nucleus to the exponential profile. 




Fig. 13. — The X-ray jet of 4C 15.55. The gray-scale ACIS image has been adaptively smoothed 
with the CIAO csmooth routine, and the green contours are 20 cm radio emission data from the 
VLA A-array observations of Hintzen, Ulvestad, & Owen (1983). The cross marking the central 
radio component has been centered on the peak of the X-ray emission. The scale can be judged 
from the pixels, which are 0'.'5 square. 

two possibilities for this excess: (1) a physical clustering of early-type galaxies around the quasars, 
and (2) massive clusters along the line of sight to the quasars, which bias the quasar catalogs via 
gravitational amplification. Our spectroscopic confirmation that at least one of the OGs in each of 
our fields has a redshift agreeing with that of the quasar to within a few hundred km s~^ supports 
the first interpretation. 

Our attempt to detect direct evidence for clusters in the two fields we discuss in this paper 
has produced ambiguous results. While we clearly seem to find an enhanced surface density of 
faint galaxies in the field of TXS 0145+386, our Chandra ACIS observation of the same field gives 
only an upper limit to any diffuse X-ray emission from the cluster. On the other hand, in the field 
of 4C 15.55, we appear to have a ~ 2.4 a X-ray detection of diffuse emission, but we cannot be 
absolutely sure that this emission is not associated with the X-ray jet from the quasar, in spite of 
our effort to excise any such emission. Neither the upper limit in the first case, nor the possible 
detection in the second, are inconsistent with the presence of a fairly massive cluster surrounding 
the respective quasars. 

Both of the OGs for which we have AO imaging appear to have compact nuclear components. 
After removal of these nuclear components, both of these galaxies are well fit by exponential profiles 
(and not by r^/'^-law profiles); however, the discrimination among profiles is much better for 4C 15.55 
ER2 because of the much higher S/N of the Subaru image. For this galaxy, there is no significant 
doubt that the stellar distribution has close to an exponential profile; even if one were to prefer the 
n = 2 Sersic model (which eliminates the need to add a compact nucleus) , it is only slightly different 
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from an exponential beyond a radius of 0'/2 (= 1.7 kpc), and, to the extent it differs, it is a worse 
fit to the data. It is true (as the referee has emphasized) that elHptical galaxies in the present-day 
universe can have a wide range of Sersic indices, including some as low as 2 (Caon, Capaccioli, &: 
D'Onofrio 1993; Khosroshahi et al. 2000). However, the same studies that have demonstrated this 
range have also shown that there is a strong positive correlation between the Sersic index and both 
the central surface brightness and the total luminosity of the galaxy. Ellipticals with n ~ 2 are, 
virtually without exception, galaxies with luminosities considerably below L* and with low central 
surface brightnesses (see, e.g.. Table 2 of Caon, Capaccioli, &: D'Onofrio 1993). There are good 
reasons to believe that they have had different formation histories than have luminous ellipticals 
(Khosroshahi et al. 2000), possibly involving secular evolution from disks. Such galaxies are in 
any case quite different from 4C 15.55 ER2, which is luminous (~ 2L*) and has a high surface 
brightness. 

Some massive galaxies strongly dominated by old stellar populations at high redshifts clearly 
show closely r^/^-law profiles (Bunker et al. 2003; Cimatti et al. 2004; Fu, Stockton, & Liu 2005; 
Daddi et al. 2005); others just as clearly show essentially pure exponential profiles (Stockton, 
Canalizo, &: Maihara 2004; Fu, Stockton, & Liu 2005; Daddi ct al. 2005). The only reasonably 
large structures in the present-day universe that show exponential profiles are the disk components 
of spiral and SO galaxies and the "pseudo-bulges" (e.g., Kormendy & Kennicutt 2004) found in 
late-type spirals. The latter are not very massive, and they are produced over a significant fraction 
of a Hubble time via secular evolution from disks. For these reasons, they are almost certainly not 
models for the exponential profiles we see in high-redshift OGs. This means that the exponential 
profiles in these massive high-redshift systems are due either to disks or to some configuration that 
is not found among massive galaxies at the present epoch. In at least one case (Stockton, Canalizo, 
&; Maihara 2004), the galaxy has a small axial ratio and looks very much like a highly inclined disk. 

If these really are massive ('^ 2 x IO^'^^Mq) disks of old stars, then they are significantly 
smaller than most massive disks at the present epoch. For example, if we take from the study of 
a complete sample of nearby disks by van der Kruit (1987) those with disk M < —20.2 (for our 
assumed Hq = 70 km s"-*^ Mpc~^), the average exponential scale length is 5.5 it 2.7 kpc, which 
corresponds to an effective radius = 9.2 ± 4.5 kpc. In comparison, TXS 0145-1-386 ERl and 
4C 15.55 ER2 have re = 1.69 kpc and 3.0 kpc, respectively. Even the larger of these is barely more 
than half as large as the smallest of the 12 luminous galaxies from the van der Kruit (1987) sample. 
Daddi et al. (2005) have similarly found that early-type galaxies at high redshifts are significantly 
smaller than their counterparts in the present-day universe. Eight out of their 9 examples were 
found to have morphologies consistent with spheroids, and only 1 appeared to be a disk of old stars. 
But it is perhaps significant that this disk-like galaxy, at z = 2.47, was by a fair margin the highest 
redshift galaxy in their sample. Daddi et al. (2005) point out that this galaxy is quite similar to 
the galaxy discussed by Stockton, Canalizo, &: Maihara (2004), which is also at a similar redshift. 

Such massive disks would carry certain implications for the early formation of massive galaxies: 
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1. Such a disk must have formed quasi- monohthically and with strong dissipation. It is unlikely 
to have suffered any major mergers after a substantial fraction of the stars had been formed. 

2. Depending on the age of the stellar population in the examples observed so far, quite large 
and sustained star formation rates are required to form most of the stellar mass in a short 
time. At a minimum, star formation rates of a few hundred Mq per year sustained for a 
few 10^ years are indicated. Note that similar constraints are indicated by the a-element 
enhancement of the stellar content of massive ellipticals at the present epoch (e.g., Worthey, 
Faber, & Gonzalez 1992; Worthey 1998). 

3. Since massive pure disks of old stars do not exist (or at least are extremely rare) at the 
present epoch, such galaxies must have a very high probability of being transformed into 
ellipticals via major mergers or into bulge-dominated SOs via less extreme encounters. This 
would not be surprising, given that our objects have been found in radio-source fields, which 
should be among the highest density regions in the early universe. At high rcdshifts, most 
of these regions will not yet have become virialized, so velocity dispersions in sub-clumps 
may well be conducive to merging activity. The fact that we suspect TXS 0145+386 ER2 
to be an ongoing merger gives some observational credence to this line of reasoning. This 
scenario might also explain why, among the admittedly small number of OGs that so far have 
published high-resolution morphological data, spheroids seem to be more common at z ~ 1.5 
than at z ~ 2.5. 

All of this fits in rather well with what we know or suspect about massive cluster ellipticals at 
the present epoch. Numerous lines of evidence (e.g., Peebles 2002; Thomas, Maraston, & Bender 
2002; Thomas et al. 2005) point to early and rapid formation of the great majority of the stars in 
the most massive and luminous ellipticals, yet some simulations (e.g., Thomas & Kauffmann 1999; 
Kauffmann & Haehnelt 2000) suggest late {z between 1 and 2) assembly of the galaxies themselves. 
The only way these two conditions can be reconciled is if the assembly is done from galaxies that 
formed most of their stars at much earlier times and that are largely free of gas. Furthermore, the 
assembly process must be dominated by (at most) a few large such galaxies in order to preserve 
the well-known color-luminosity relation for ellipticals. Thus, early forming massive disks in which 
star formation has been driven to completion or truncated by the expulsion of remaining gas could 
plausibly act as reservoirs for the stars that eventually end up in present-day elliptical or cD 
galaxies. 
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the presentation and content of the paper. This research has been partially supported by NSF grant 
AST03-07335 and by the National Aeronautics and Space Administration through Chandra Award 
Number GO2-3180X issued by the Chandra X-ray Observatory Center, which is operated by the 
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ministration under contract NAS8-03060. It made use of the NASA/IPAC Extragalactic Database 



- 22 - 



(NED) which is operated by the Jet Propulsion Laboratory, Cahfornia Institute of Technology, un- 
der contract with the National Aeronautics and Space Administration. The authors recognize the 
very significant cultural role that the summit of Mauna Kea has within the indigenous Hawaiian 
community and are grateful to have had the opportunity to conduct observations from it. 

REFERENCES 

Barr, J. M., Bremer, M. N., Baker, J. C, & Lehnert, M. D. 2003, MNRAS, 346, 229 
Best, P. N. 2000, MNRAS, 317, 720 

Best, P. N., Lehnert, M. D., Miley, G. K., & Rottgering, H. J. A. 2003, MNRAS, 343, 1 
Brown, M. J. L, Jannuzi, B. T., Dey, A., & Tiede, G. P. 2005, ApJ, 621, 41 
Bruzual, G. & Chariot, S. 2003 MNRAS, 344, 1000 

Bunker, A., et al. 2003, in The Mass of Galaxies at Low and High Redshift, (ESO, Garching), 
p. 262 

Calzetti, D., Armus, L., Bohlin, R. C., Kinney, A. L., Kornneef, J., & Storchi-Bergmann, T. 2000, 
ApJ, 533, 682 

Caon, N., Capaccioh, M., & D'Onofrio, M. 1993, MNRAS, 265, 1013 
Cimatti, A., et al. 2002, A&A, 381, L68 
Cimatti, A., et al. 2003, A&A, 412, LI 
Cimatti, A., et al. 2004, Nature, 430, 184 

Chapman, S. C, McCarthy, P., k Persson, S. E. 2000, AJ, 120, 1612 
Chen, H.-W., & Marzke, R. 2004, ApJ, 615, 603 

Croom, S. M., Shanks, T., Boyle, B. J., Smith, R. J. Miller, L., Loaring, N. S., & Hoyle, F. 2001, 
MNRAS, 325, 483 

Daddi, E., Cimatti, A., Pozzetti, L. Hoekstra, H., Rottgering, H. J. A., Renzini, A., Zamorani, G., 
& Mannucci, F. 2000, A&A, 361, 535. 

Daddi, E., Cimatti, A., Renzini, A., Fontana, A., Mignoli, Pozzetti, L., Tozzi. P., &; Zamorani, G. 
2004, ApJ, 617, 746 

Daddi, E., et al. 2005 ApJ, 626, 680 

Dey, A., Graham, J. R., Ivison, R. J., Small, L Wright, G. S., & Liu, M. C. 1999, ApJ, 519, 610 



-23- 



Dunlop, J., Peacock, J., Spinrad, H., Dey, A., Jimenez, R., Stern, D., & Windhorst, R. 1996, 
Nature, 381, 581 

Ebeling, H. 2003, MNRAS, 340, 1269 

Ferrara, A., Bianchi, S., Cimatti, A., &: Giovanardi, C. 1999, ApJS, 123, 437 
Fu, H., Stockton, A., & Liu, M. 2005, ApJ, 632, 831 
Gehrels, N. 1986, ApJ, 303, 336 

Gilbank, D. G., Smail, I., Ivison, R. J., k Packham, C. 2003, MNRAS, 346, 1125 
Hall, P., et al. 2001, AJ, 121, 1840 

Hintzen, P., Ulvestad, J., & Owen, F. 1983, A J, 88, 709 
lye, M. et al. 2003, ApJ, 590, 770 

Kauffmann, G., k Haelinelt, M. 2000, MNRAS, 311, 576 

Kewley, L. J., Geller, M. J., & Jansen, R. A. 2004, AJ, 127, 2002 

Khosroshahi, H. G., Wadadekar, Y., Kembhavi, A., & Mobasher, B. 2000, ApJ, L103 

Kinney, A. L., Calzctti, D., Bohlin, R. C., McQuade, K., Storchi-Bergmann, T., k Schmitt, H. R. 
1996, ApJ, 467, 38 

Kobayashi, N., et al. 2000, in Proc. SPIE 4008: Optical and IR Telescope Instrumentation and 
Detectors, ed. M. lye k A. F. Moorwood, p. 1056 

Kormendy, J., k Kennicutt, R. C., Jr. 2004, ARA&A, 42, 603 

Liu, M. C., Dey, A., Graham, J. R., Bundy, K. A., Steidel, C. C., Adelberger, K., k Dickinson, M. 
E. 2000, AJ, 119, 2556 

McCarthy, P. J., et al. 2001, ApJ, 560, L131 

Moriando, G., Cimatti, A., k Daddi, E. 2000, A&A, 364, 26 

Mulhs, C. R., Rosati, P., Lamer, G, Bohringer, H., Schwope, A., Schuecker, P., k Fassbender, R. 
2005, ApJ, 623, L85 

Oke, J. B., et al. 1995, PASP, 107, 375 

Peebles, P. J. E. 2002, ASP Conf. Ser., 283, 351 

Peng, C. Y., Ho, L. C, Impey, C. D., k Rix, H.-W. 2002, AJ, 124, 266 
Pozzetti, L., k Mannucci, F. 2000, MNRAS, 317, L17 



-24- 



Rottgering, H., De Breuck, C, Daddi, E., Kurk, J., Miley, G., Pentericci, L., Overzier, R., & 
Venemans, B. 2004, 

Saracco, P., Longhetti, M., Giallongo, E., Arnouts, S., Cristiani, S., D'Odorico, S., Fontana, A., 
Nonino, M., & Vanzella, E. 2004, A&A, 420, 125 

Sawicki, M., Stevenson, M., Barrientos, L. F., Gladman, B., Mallen-Ornelas, G., & van den Bergh, 
S. 2005, ApJ, 627, 621 

Smith, G. P., Smail, I., Kneib, J.-P., Davis, C. J., Takamiya, M., Ebeling, H., & Czoske, O. 2002b, 
MNRAS, 333, L16 

Smith, H. E., Burbidgc, E. M., Baldwin, J. A., Tohhne, J. E., Wampler, E. J., Hazard, C, k 
Murdoch, H. S. 1978, ApJ, 215, 427 

Soifer, B. T., Matthews, K., Neugebauer, G., Armus, L., Cohen, J. G., Persson, S. E., k, Smail, I. 
1999, AJ, 118, 2065 

Spinrad, H., Dey, A., Stern, D., Dunlop, J, Peacock, J., Jimenez, R., k Windhorst, R. 1997, ApJ, 
484, 581 

Stiavelh, M., k Treu, T. 2001, in Galaxy Disks and Disk Galaxies, ASP Conf. Ser., 230, ed. Funes 
k Corsini, (ASP, San Francisco), p. 603 

Stockton, A. 2001, in Astrophysical Ages and Timescales, ASP Conf. Ser., 245, ed. von Hippel et 
al., (ASP, San Francisco), p. 517 

Stockton, A., Canahzo, G., k Maihara, T. 2004, ApJ, 605, 37 

Stockton, A., Kellogg, M., k Ridgway, S.E. 1995, ApJ, 443, L69 

Thomas, D., k KaufFmann, G. 1999, in Spectrophotometric Dating of Stars and Galaxies, ASP 
Conf. Ser., 192, ed. I. Hubeny, S. R. Heap, k R. H. Cornett, p. 261. 

Thomas, D., Maraston, C, k Bender, R. 2002, in Reviews in Modern Astronomy, ed. R. E. 
Schielicke, v. 15, p. 219 [astro-ph/0202166] 

Thomas, D., Maraston, C, Bender, R., k de Oliveira, C. M. 2005, ApJ, 621, 673 

van der Kruit, P. C. 1987, A&A, 173, 59 

van Dokkum, P. G., k Ellis, R. S. 2003, ApJ, 592, L53 

Wills, D., k Lynds, C. R. 1978, ApJS, 36, 317 

Wold, M., Armus, L., Neugebauer, G., Jarrett, T. H., k Lehnert, M. D. 2003, AJ, 126, 1776 
Worthey, G. 1998, PASP, 110, 888 



-25- 



Worthey, G., Faber, S. M., & Gonzalez, J. J. 1992, ApJ, 398, 69 
Yan, L. et al. 2004b, ApJS, 154, 75 

Yan, L., Thompson, D., & Soifer, B. T. 2004, AJ, 127, 1274 



This preprint was prepared with the AAS IM^^X macros v5.2. 



